Scattering of Magnetic Solitons in two dimensions 
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Solitons which have the form of a vortex-antivortex pair have recently been found in the Landau- 
Lifshitz equation which is the standard model for the ferromagnet. We simulate numerically head-on 
collisions of two vortex-antivortex pairs and observe a right angle scattering pattern. We offer a 
resolution of this nontrivial dynamical behaviour by examining the Hamiltonian structure of the 
model, specifically the linear momentum of the two solitons. We further investigate the dynamics 
of vortices in a modified nonlinear cr-model which arises in the description of antiferromagnets. We 
confirm numerically that a robust feature of the dynamics is the right angle scattering of two vortices 
which collide head-on. A generalization of our theory is given for this model which offers arguments 
towards an understanding of the observed dynamical behaviour. 
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I. INTRODUCTION 

Localized solutions, often called solitons, play an in- 
creasingly important role in nonlinear field theories in 
two dimensions. Topological structures exist in partic- 
ular in magnetic systems and have been studied exten- 
sively, both theoretically and experimentally 

An easy-plane ferromagnet is described by the Landau- 
Lifshitz equation 



h = n X f , 



(1) 



f = An 



nse. 



1. 



The field n represents the local magnetization of the ma- 
terial. We shall study here the case of a two-dimensional 
medium so we assume n = n(a;,y,t). The dot denotes a 
time derivative, rt3 is the third component of n, A is the 
Laplace operator and e — (0, 0, 1) is the unit vector in 
the third direction. The normalization condition n = 1 
imposed in the initial condition is preserved by the equa- 
tions of motion. 

Static solutions of model (|l|) are the well-studied vor- 
tices. Isolated vortices are spontaneously pinned objects 
that is no vortex in free translational motion can be found 
in a 2D ferromagnet. The same is true for any isolated 
object with nontrivial topology in a 2D ferromagnet 
the most well-known example being the magnetic bubbles 
in easy-axis ferromagnetic films 

Coherently traveling solutions of (0) have been found 
in ||]. They have the form of a vortex-antivortex pair 
and their velocity may take any value between zero and 
unity, which is the velocity of magnons in the system. 

We now turn to a different class of systems, namely 
antiferromagnets. The dynamics of the staggered mag- 
netization in the antiferromagnetic continuum is given by 
the nonlinear cr-model 



f = An — e. 



n^ = l. 



n X [ii — f] =0, 



(2) 



where the double dot denotes a second time derivative. 
The above model has the same static vortex solutions as 
(|l|). On the other hand, vortices in model (^) can be 
found in free translational motion. This is due to the 
fact that the model is invariant under Lorentz transfor- 
mations. 

Our main purpose is to study collisions of solitons in 
models (|l|) and (g). In the ferromagnet no collision be- 
tween two vortices can take place. Two vortices with the 
same topological charge will rotate around each other 
while a vortex and an antivortex undergo Kelvin mo- 
tion perpendicular to the line connecting them. However, 
collisions can occur between two vortex-antivortex pairs. 
We elaborate on the arguments of and argue that 
head-on collisions between vortex-antivortex pair solitons 
give a right angle scattering pattern. 

We also study collisions between vortices in antiferro- 
magnets in Eq. (|^). They scatter at right angles as found 
in numerical simulations. In fact, the right angle scatter- 
ing phenomenon seems to be a robust feature in vari- 
ous two-dimensional models which have soliton solutions 
1^,^. However, it is a nontrivial and strange behaviour 
at least from the point of view of scattering of ordinary 
particles. 

The colliding objects in the two models that we study 
are essentially different from each other. While vortices 
within model (H) are topologically nontrivial objects, the 
colliding vortex-antivortex pairs in have a vanishing 
topological charge. However, we argue that the under- 
lying Hamiltonian structure allows to study the soliton 
interaction in the two models in close analogy. 

The outline of the paper is as follows. In Section II we 
simulate head-on collisions of vortex-antivortex pairs in 
(|l|) and give a theoretical description which exploits the 
form of the linear momentum. In Section III the results 
of head-on collision simulations of vortices in model (|^) 
are given together with arguments for the understanding 
of this behaviour. The conclusions are given in Section 
IV. 
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II. HEAD-ON COLLISIONS OF 
VORTEX-ANTIVORTEX PAIRS IN PLANAR 
FERROMAGNETS 



A ferromagnet can be described in terms of a magneti- 
zation vector which satisfies the Landau-Lifshitz equation 
(H) . The constraint on the field n can be resolved and the 
theory can be formulated in terms of a complex variable 



1 + n3 



which satisfies the equation 

i VL = -Arj + = d^n dM = n. 



(3) 



(4) 



O denotes the complex conjugate of 0. 

We use the formulation through the complex variable 
O in all numerical simulations. We avoid the formula- 
tion through the vector variable n since the constraint 
on it makes an accurate computer calculation of the time 
derivatives of the field rather cumbersome. 

The model has some interesting static vortex solutions 
of the form 



r!° = /(p) e-^ 



K = ±1, 



(5) 



where p, 4> are polar coordinates and f{p = 0) = 0, f{p 
oo) 1. We call the configuration with k = 1 a vortex 
and the one with k — —1 an antivortex. Vortex solutions 
have infinite energy and it has been argued that they are 
physically relevant [^|6| . 

In the study of the dynamics of magnetic vortices the 
central roleJ_s played by a scalar quantity called the local 
vorticity 



7 = d^TT (9^V, 



(6) 



where e^^ is the two-dimensional totally antisymmetric 
tensor. The two components of the linear momentum are 
then expressed as 



y-f dxdy, 



Py — 1^1 dxdy. 



(7) 



Of fundamental importance is the Poisson bracket re- 
lation between the two components of the linear momen- 
tum. This reads 



{Px,Py} 



(8) 
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FIG. 1. A simple vorticity distribution of a soliton con- 
sists of two lumps with opposite sign. They are here symmet- 
rically placed on either side of the a;-axis. The lower shaded 
area represents positive vorticity while the upper shaded area 
represents negative vorticity. P denotes the linear momentum 
of the pair. 



where 



T — J dxdy 



(9) 



is the total vorticity. 

In the present model n — cosQ and ^ = <I> have been 
used Q as the canonical fields. They arc defined through 
rii ~ cos8 sin$, n2 — cosQ sin$, 713 = sin8. The explicit 
form of the vorticity is 



7 = e^i/ sin6 9^9 d^^. 
The total vorticity of a vortex is 



r 



7 dxdy = —2 TTK, 



(10) 



(11) 



that is, r = — 2 7r for vortices and r = 2 7r for antivortices. 

The implications of a nonvanishing total vorticity to 
the dynamics is an issue which has been thoroughly stud- 
ied in the case of magnetic vortices and bubbles [p|,p|JTl|] as 
well as in the case of vortices in other interesting mod- 
els p^ , p^ . The most striking result is that it leads to 
spontaneous pinning of these topological objects. 

Since a nonvanishing total vorticity T implies pinning 
of an object, we infer that a solution which moves freely 
should have a vanishing F. In this respect, the vortex- 
antivortex ansatz offers the simplest possibility. Fig. |l| 
gives a schematic representation of it. This consists of 
two lumps, one having negative and the other one posi- 
tive sign. We suppose that the vortex is roughly laying in 
the shaded area with the negative sign and the antivor- 
tex in the shaded area with the positive sign. This figure 



2 



is supposed to act only as a guide for our discussion and 
there is no strict way to distinguish the two vortices and 
define their positions. However, if relation (^) is applied 
to each vortex separately then the quantity on the right 
hand side is nonvanishing. It is then implied that each 
vortex will propagate in the horizontal direction under 
the influence of the other vortex. The picture is consis- 
tent with linear momentum considerations. That is, an 
application of Eq. to the full ansatz gives a nonvan- 
ishing x-component of the linear momentum. Fig. |l| will 
serve in the following discussion as a prototype and will 
motivate our theoretical arguments. 

The Kelvin motion of a vortex-antivortex pair in a fer- 
romagnet has been investigated in [ p^ . The motion of a 
bubble- antibubble ansatz has also been studied |l^. The 
situation is found to be similar in some other systems 
such as an antiferromagnet immersed in a uniform mag- 
netic field 1^, a model for superconductors jl^, for su- 
perfluid helium [ pTf and the nonlinear Schrodinger equa- 
tion . It has been pointed out that the gross features 
of this dynamical behaviour are analogous to the planar 
motion of charges under the influence of a magnetic field 
perpendicular to the plane. In particular, two oppositely 
charged particles undergo Kelvin motion traveling along 
parallel trajectories. The analogy has been made precise 
by use of relation (|^) and an analogous relation in the 
charge motion problem 

The calculation of steadily moving coherent structures 
in a 2D ferromagnet, which have the form of a vortex- 
antivortex pair has been done in We have used here 
the numerical code of to reproduce them since there 
is no available analytical formula. Fig. ^ is an example 
contour plot for a soliton with velocity v = 0.5. The up- 
per entry gives contour plots of the quantity 10 |ri| and 
the two-vortex character of the configuration is rather 
obvious. The lower entry is a contour plot for the local 
vorticity ( |l0[) . An important result of the analysis in |^ is 
that the velocity is coUinear with the linear momentum. 

We are now sufficiently motivated to explore the possi- 
bility of scattering of vortex-antivortex pair solitons. We 
denote by A„(a;, y) the solution with velocity v along the 
X-axis (set t = 0). The product ansatz 

n{x,y)^A,(^x+^,y^ A^^,(^x-^,y^ (12) 

represents two vortex-antivortex pair solitons at a dis- 
tance S apart which are in a head-on collision course. The 
ansatz ( |l2| ) is used as an initial condition in a straight- 
forward numerical integration of Eq. (^. We typically 
set v = 0.5, 5=10. 

We have set up a numerical mesh as large as 600 x 600 
with uniform lattice spacing h = 0.1. The space deriva- 
tives are calculated by finite differences and the time in- 
tegration is performed by a fourth order Runge-Kutta 
method. The results are presented in Figs. ^ and^. Fig. ^ 
presents a contour plot for the field 10 at three char- 
acteristic snapshots. In the first entry the initial ansatz 



( [121 ) is shown. The second snapshot is taken when the 
solitons are more or less at a minimum separation. No 
vortex-antivortex annihilation process takes place. This 
behaviour should be expected since the vortex-antivortex 
pair solitons are stable solutions of the equation. The ar- 
gument is supported by numerical simulations showing 
that a vortex-antivortex ansatz preserves its character 
when traveling, provided that the vortex and antivortex 
are not very close to each other 0,^. The last snap- 
shot shows the system after the collision. A right angle 
scattering pattern has been produced. 
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FIG. 2. Contour plot for a vortex-antivortex pair soliton 
in a ferromagnet with velocity v = Q.5. The upper entry gives 
contour plots of the quantity 10 |f2|. We plot the levels 1, 3, 
5, 7, 9. The lower entry is a contour plot of the local vorticity 
for the same soliton. Solid lines represent positive values and 
dashed-dotted lines negative values of vorticity. We plot the 
levels ±0.1, ±0.2, ±0.4, ±0.8, ±1.2. 
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FIG. 3. Contour plot of the field 10 \ at three charac- 
teristic snapshots of the head-on collision simulation of vor- 
tex-antivortex pair solitons in ferromagnets. It is shown: the 
initial ansatz (upper entry, time t — Q), a snapshot at the time 
of collision (middle entry, t = 6.6), and well after collision 
(lower entry, t — 13.2). Contour levels as in Fig. |^. 



FIG. 4. Contour plot of the local vorticity 7 of Eq. 
for the solitons of Fig. M. Contour levels as in Fig. ti. 
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FIG. 5. Stars denote the zeros of Q during the numerical 
simulation of Fig. ^ We also trace the maximums (circles) 
and minimums (diamonds) of the vorticity. The solitons are 
initially located at AB and CD respectively. Symbols are 
plotted every 0.4 time units. 



The situation becomes clearer in Fig ^ where we rep- 
resent the solitons in terms of their local vorticity dis- 
tribution. The soliton on the left half plane should be 
compared directly with that in the lower entry of Fig. |^. 
It obviously has a linear momentum and velocity point- 
ing to the positive x-direction. The other soliton has the 
opposite linear momentum and velocity. 

It is rather clear from the picture that the solitons will 
not bounce back after collision. This is precluded by the 
form of the local vorticity distribution. This possibility 
would require that the vortex and antivortex interchange 
their position. Instead, the possibility appears that, at 
collision time, the two pairs of vorticity lumps in the 
upper and lower half-planes will form two new vortex- 
antivortex pairs. 

One has to apply Eq. for each of the vorticity lumps 
separately. Alternatively, one can consider pairs of lumps 
which tend to travel parallel to each other undergoing 
Kelvin motion. Application of this idea to Fig. deter- 
mine the time evolution of the system. Finally, the two 
pairs on the upper and lower half planes tend to travel 
parallel to each other along the y-axis and form bound 
states. 

An equivalent point of view is to follow the linear mo- 
mentum of each soliton separately. The linear momen- 
tums of the outgoing solitons clearly lay on the y-axis 
and have opposite sign. 

A subtle but important question is whether we can ap- 
ply Eq. (H) separately for each of the vortices which con- 
sist the vortex-antivortex pair. A rigorous answer can 
not be given here. On the other hand the construction of 
solitary waves in [p|j7||r^ suggests an affirmative answer 



whose range of validity is interesting to study. 

The two solitons emerging after collision are very sim- 
ilar to the initial ones though not exactly the same. In 
fact the drift velocity of the outgoing solitons is some- 
what larger. In Fig. || we have traced, during the time 
evolution, the points where the complex field ft vanishes 
and also the points where the vorticity 7 attains its max- 
imum and minimum values. The two kinds of extrema 
are close during the whole period of time evolution. This 
is because the solitons used in the simulations of this 
chapter have a pronounced vortex-antivortex character. 
In |1^ traveling solutions of the Landau-Lifshitz equa- 
tion have been studied which are different than the ones 
used here. Numerical simulations of scattering of these 
solitons also produce a right angle pattern. 

It is possible to give a picture, corresponding to the 
scattering of vortex-antivortex pairs, in terms of 2D mo- 
tion of charged particles interacting via their electric field 
and placed in a magnetic field perpendicular to the plane. 
In fact, we have to consider two electron-positron pairs. 
Consider the first electron-positron pair located at points 
A, B of Fig. ||. and the second pair at C, D. The charges 
move similar to the vortices. Their actual orbits resemble 
those for the solitons shown in Fig. ||. 

Our last remark in this section goes to some related 
work in hydrodynamics. There are solutions of the 
two-dimensional Euler equations which describe vortic- 
ity dipoles. Note that, in this context, vorticity has its 
ordinary hydrodynamic meaning. The best-known such 
solution seems to be the Lamb dipole [^. Another one 
has been found in |^ . A head-on collision between two 
dipoles produces a pattern analogous to that in Fig. ^ 
of the present paper |pi|,P^ . Furthermore, a simple con- 
struction is given in [pO[ page 223, to which our Fig. || 
can be compared. Further interesting cases of scatter- 
ing between pairs of objects in hydrodynamics have been 
studied. The most complex behaviour has been observed 
in 1 23 1 and includes stochastic and quasiperiodic motion 
of vortices. 



III. HEAD-ON COLLISIONS OF VORTICES IN 
ANTIFERROMAGNETS 

Our main objective is to study scattering of vortices 
within model and to show that the process can 
be studied in close analogy to the corresponding phe- 
nomenon in the ferromagnet. A right angle scattering 
behaviour of solitons has been observed in the isotropic 
(7- model 1^, that is model (||) without the anisotropy 
term. 

The examination of the local vorticity 7 has led to a 
successful approach for the collision of vortex-antivortex 
pairs in ferromagnets in Section II. We find it instructive 
to look at the collision process in terms of the vorticity 
also in the present model. A simple generalization of def- 
inition can be used . The vorticity attains a simple 
form when it is expressed in terms of the vector field n: 
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7 = df,h ■ d^n = e^^, 5^(h • d^n). (13) 

In the equation for an antiferromagnet in a uniform 
magnetic field was studied. Vortices in this system are 
spontaneously pinned, thus their dynamics is analogous 
to that of ferromagnetic vortices. This unexpected be- 
haviour is probed by a topological term which enters the 
vorticity. However, such a term is absent in the model 
studied in this section. 

The vorticity (^3|) has the form of a total divergence 
and can be integrated in all space to show that the total 
vorticity vanishes for solutions with reasonable behaviour 
at infinity: 



7 dxdy — 0. 



(14) 



In particular, it vanishes for the vortex solutions. Re- 
lations (^ - (||) apply in the present context without 
modification and they will be the fundamental relations 
to be used in the following analysis. 

Eq. (13) shows that for a static vortex 7 vanishes iden- 
tically. On the other hand, we can obtain a steadily 
traveling vortex by applying a Lorentz transformation to 
the static vortex solution (||). We denote the traveling 
vortex by ^1° and the velocity is < u < 1. The distri- 
bution of 7 for a Lorentz boosted vortex is nonvanishing 
and can be calculated numerically. The vortex with ve- 
locity v = 0.7 is represented by a contour plot of the field 
10|f2| in the upper entry of Fig. ^. A corresponding plot 
for 7 is given in the lower entry of the figure. 

The vorticity distribution has the form of two lumps, 
thus it resembles the sketch of Fig. |l|. This is no sur- 
prise. In fact the following two remarks make it plausible. 
Firstly, we see that the total vorticity vanish according 
to relation (|lj). Secondly, an inspection of the form (|7) 
of the linear momentum makes it clear that a nonvan- 
ishing component is furnished by two lumps of vorticity 
with opposite signs. This is not the only form of local 
vorticity that furnishes a nonvanishing linear momentum 
but it is certainly the simplest. Since the vortex solution 
with K = 1 is indeed the one with the simplest topological 
complexity, we expect its vorticity distribution to have 
the simplest possible form. 

We calculate numerically the points where the maxi- 
mum and minimum of the vorticity lumps are located. 
It turns out that these points are the (0, ±0.59) for any 
value of the velocity v. 

A further example on the present ideas is offered by 
the Belavin-Polyakov solutions . We apply a Lorentz 
transformation, with velocity u, to the simplest one: 

n 



(x — vt) I \J\ — + iy. Its local vorticity is 
16 V 



7 : 



V 



1 



1 



(X 



Vt)' 



v 



(15) 



In accordance with the above remarks, it has the shape 
of two lumps with opposite sign, located on either side 
of the X-axis and traveling along the a;-axis. 
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FIG. 6. Contour plot for a traveling vortex with velocity 
V — 0.7. The upper entry gives contours of the field 10 |n|. 
The lower entry is a contour plot of the local vorticity 7. 
Solid lines represent positive values and dashed-dotted lines 
represent negative values of vorticity. Contour levels as in 
Fig. I 



We are now ready to present numerical simulations of 
head-on collisions of vortices. We make an ansatz repre- 
senting two vortices. The choice is not unique and the 
simplest one seems to be the product ansatz: 



n{x,y)^n° (x+-,y] n" 



(16) 



where $7° is the single vortex solution given in Eq. (||). 
The two vortices are a distance S apart. 

The numerical mesh as well as the details of the algo- 
rithm that we use here are similar to those of Section II. 
We use vortices with k = 1. They are initially at rest but 
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immediately start to drift away from each other due to 
their mutual repulsion and escape to infinity. 

In order to invoke a head-on collision between vortices 
we consider the product ansatz of two vortices which have 
opposite velocities: 



(17) 



i^y{x, y) denotes the Lorentz transformed vortex solution 
with velocity -y, at t — Q. 

We typically use (5 = 6. In all simulations the vortices 
start to move against each other with velocities close to 
the value v but they immediately begin to decelerate due 
to their mutual repulsion. The future of the process de- 
pends crucially on the magnitude of the velocity. At low 
velocities the two vortices approach to a minimum dis- 
tance at which they come to rest and then turn round 
and move off in opposite directions. When the velocity 
exceeds a critical value (which is i;c~0.65 for 5 = 6) the 
two vortices collide and scatter at right angles. This re- 
sult does not depend on the details of the initial ansatz or 
on the initial velocity of the vortices, as long as this ex- 
ceeds the critical value. We have tested our algorithm for 
velocities up to the value w = 0.9 in ansatz (p^. However, 
one must keep in mind that the velocity of the vortices 
at the time of collision is smaller than the velocity in the 
initial ansatz. 

In Fig. 1^ we give a contour plot for the norm of the 
field at three characteristic snapshots. The first en- 
try presents the initial configuration (p7|). In the middle 
snapshot, taken at collision time, it is clear that the two 
vortices come on top of each other. There is no topologi- 
cal reason, related to the field fi, that could prevent this 
double vortex to form and there is also no such reason 
that could prevent the vortices either to continue trav- 
eling in the horizontal direction or to reemerge traveling 
in the vertical direction. We add that, at the present 
level of description, we can find no reason that would 
enforce them to follow either of the two possibilities. In 
the last snapshot the two new vortices that emerge after 
the collision, are drifting away from each other along the 
y-axis. 

We find it instructive to look at the collision process 
using the vorticity. Our description will closely follow 
that in Section II in connection with the scattering of 
vortex-antivortex pairs. The dynamics in both systems 
is determined by the corresponding vorticity distribution. 
A comparison of the lower entries of Figs. p| and |^ gives a 
hint that the underlying dynamics should be of a similar 
nature in both models. 

In Fig. 1^ we give the vorticity at three snapshots which 
correspond to those of Fig. 0. Fig. || should be compared 
directly with Fig. ^. An examination of these results 
shows that the arguments of Section II for the soliton 
scattering which rely upon the linear momentum rela- 
tions (0), (^) are apphcable here, too. 
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FIG. 7. Contour plot of the field 10 |f2| at three character- 
istic snapshots for the head-on collision simulation of vortices 
in an antiferromagnet. It is shown: the initial ansatz (upper 
entry, time t = 0), a snapshot at the time of collision (middle 
entry, t = 4.2), and well after collision (lower entry, t = 8.4). 
Contour levels as in Fig. B. 
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FIG. 8. Contour plot of the local vorticity 7 of Eq. for 
the vortices of Fig. ^ Contour levels as in Fig. ^. 

The upper entry in the figure corresponds to the ini- 
tial ansatz and each of the two vorticity dipoles should be 
compared to that given in the lower entry of Fig. ^. The 
two vortices in the first entry of Fig. |^ are approaching 
each other while their dynamical features, as described 



by the vorticity, are not substantially modified. The re- 
pulsion which could decelerate them and make them turn 
round, is overcompensated by the large enough initial ve- 
locity. When the two vortices come close to each other 
(middle entry) the two pairs of vorticity lumps, lying in 
the upper and lower half plane, interact. The subsequent 
evolution of the vorticity lumps is governed by relation 
(^. In particular, this relation has to be applied to each 
of the vorticity lumps separately. As is indicated by the 
simulation, the lumps survive throughout the process and 
the simple dynamics implied by is sustained during 
the scattering process. Following the discussion in Sec- 
tion II we examine the dynamics of lumps in a pairwise 
manner. The two pairs in the upper and lower plane 
tend to move along the y-axis. Consequently, the initial 
partners separate and two new vortices are formed which 
travel in opposite directions along the y-axis, as shown 
in the lower entry of the figure. 

An equivalent approach is to follow the linear momen- 
tum of the solitons. The two pairs in the lower and upper 
half-plane have their linear momentums lying along the 
y-axis but with opposite signs. They subsequently tend 
to go off along the y-axis. The important numerical re- 
sult is that the vorticity lumps roughly preserve their 
shape throughout the process. This is due to the fact 
that traveling vortices are stable solutions of the model. 

We note here that our result is obtained when Eq. (||) 
is applied to each vorticity lump separately. We have 
been motivated to follow this approach because of its 
success with respect to studying the dynamics of vortex- 
antivortex pairs in a ferromagnet and because of the con- 
sistency of the picture with the linear momentum con- 
siderations. On the other hand a rigorous proof of its 
validity is lacking. 
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FIG. 9. Stars denote the zeros of Q during the numerical 
simulation of Fig. |^. We also trace the maximums (circles) 
and minimums (diamonds) of the vorticity. The vortices are 
initially located at points A and B. Symbols are plotted every 
0.4 time units. 
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In Fig. [q] we track some characteristic points of the 
vortices throughout the simulation time. The stars de- 
note successive points where the centers of the two vor- 
tices he during the simulation, that is, where the complex 
function vanishes. The circles denote successive loca- 
tions of the maximum and the diamonds the locations 
of the minimum of the vorticity distribution. We plot 
two circles and two diamonds at every time instant. At 
the beginning of the simulation the two vortices are cen- 
tered at points A and B respectively. They immediately 
start to decelerate but when their centers reach a dis- 
tance « 2.4 they seem to accelerate considerably, they 
eventually merge at the origin and then they separate 
along the y-axis. The trajectories of the extrema of the 
vorticity show that after the collision it takes some time 
until the two new vortices are organized again. The ve- 
locity at the end of the numerical simulation is somewhat 
lower than that in the initial ansatz. This should be due 
to spin waves emitted during the process. 

The remarks of the last paragraph on the motion of the 
vortex centers are in agreement with an analytical solu- 
tion obtained in |25[ | representing scattering of solitons in 
an integrable chiral model. Close to collision time, the 
centers move according to the law x « ±a/^ {t < 0). 
This gives a velocity dx/dt ^ oo as t ^ Q. They collide 
at t = and the centers of the two new solitons, which 
emerge along the y axis, follow y w ±^/^ (t > 0). 

The arguments presented in this section are certainly 
not sufhcient to exclude, other than right angle, interest- 
ing possibilities of interaction and scattering of solitons. 
Nevertheless, they imply that the right angle scattering 
process is expected to be generic for solitons in two- 
dimensional Hamiltonian models. In the present work we 
have found no relation of the topology of the solitons to 
the right-angle scattering phenomenon, Therefore right- 
angle scattering is also expected to occur among non- 
topological solitons [E6|. 



IV. CONCLUSIONS 

We have given a description of the right angle phe- 
nomenon of solitons through numerical simulations, as 
well as arguments which suggest that it should be generic 
in two dimensions. Two systems have been examined. 
Vortex-antivortcx pairs in planar ferromagnets and vor- 
tices in antiferromagnets. The peculiar scattering be- 
haviour is mainly attributed to the fact that solitons are 
extended structures rather than point like particles. This 
point is accounted for by the representation of a travel- 
ing soliton through a pair of lumps. Furthermore, we find 
that these two lumps act as independent physical entities 
at the time of collision. 

It is desirable to observe experimentally scattering of 
solitons in ferromagnetic and antiferromagnetic films. 
We expect that the present theoretical analysis will be 
useful in studies of systems of a lot of vortices and in 



particular in studies of the thermodynamics of magnetic 
systems. Suffice it to say that, in a magnetic material, 
vortices are expected to appear in pairs. In the study of 
the thermodynamics of layered antiferromagnets one ex- 
pects to find the signature of topological excitations. The 
remark may prove important especially in view of the dif- 
ficulty to observe isolated antiferromagnetic solitons due 
to the lack of a significant total net magnetization. 

The right angle scattering pattern appears also and 
has been understood from the point of view of the geom- 
etry of the moduli space for BPS monopoles in a three- 
dimensional Yang-Mills-Higgs theory j27j . 

We have proceeded in the simulation of scattering of a 
vortex and an antivortex in the cr-model . We use an 
initial ansatz similar to that of Eq. ( p^ ) . We simulate the 
time evolution of the system numerically and find that 
the vortices attract each other. They eventually collide 
at the origin and annihilate. It is quite interesting that 
the energy is dissipated at right angles. The phenomenon 
is presumably closely related to the present study. A sim- 
ilar simulation with corresponding results has been done 
in 1 28 1 for non-gauged cosmic strings. 

The basic dynamical structure which leads to right an- 
gle scattering in Hamiltonian models is also present in 
the complex Ginzburg-Landau equation (CGLE) which 
describes nonlinear oscillatory media . There is actu- 
ally a variety of nonconservative systems where scattering 
behaviour analogous to that studied in the present paper 
has been observed. An interesting example is a 2D fluid 
layer subjected to externally imposed oscillations. Co- 
herent structures are formed which are experimentally 
observed to scatter at an almost right angle M . 
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